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^| Observations of the thermal X-ray emission from radio pulsars implicate that the size of hot spots 

(S| is much smaller then the size of the polar cap that follows from the purely dipolar geometry 

of pulsar magnetic field. Most plausible explanation of this phenomena is an assumption that 
ry^ the magnetic field at the stellar surface differs essentially from the purely dipolar field. We can 

jT^ determine magnetic field at the surface by the conservation of the magnetic flux through the 

area bounded by open magnetic field lines. Then the value of the surface magnetic field can be 
estimated as of the order of 10 14 G. On the other hand observations show that the temperature of 
I the hot spot is about a few million Kelvins. Based on these observations the Partially Screened 

^ Gap (PSG) model was proposed which assumes that the temperature of the actual polar cap (hot 

spot) equals to the so called critical temperature. 
, We discuss correlation between the temperature and corresponding area of the thermal X-ray 

emission for a number of pulsars. The results of our analysis show that the PSG model is suitable 
to explain both cases: when the hot spot is smaller and larger then conventional polar cap. We 
argue that in the second case structure and curvature of field lines allow pair creation in the closed 
\Q field lines region thus the secondary particles can heat the stellar surface outside the actual polar 

cap. 

We have found that the Curvature Radiation (CR) plays dominant role in avalanche pair produc- 

00 

tion in the PSG. We studied dependence of the PSG parameters on the pulsar period, the magnetic 
field strength and the curvature of field lines. 
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1. Introduction 

The Standard model of radio pulsars assumes that there exists the Inner Acceleration Region 
(IAR) above the polar cap where the electric field has a component along the opened magnetic field 
lines. In this region particles (electrons and positrons) are accelerated in both directions: outward 
and toward the stellar surface. Consequently, outflowing particles are responsible for generation 
of the magnetospheric emission (radio and high-frequency) while the backflowing particles heat 
the surface and provide required energy for the thermal emission. In such scenario analysis of X- 
ray radiation is an excellent method to get insight into the most intriguing region of the neutron star. 

1.1 Observations 

About 30 years ago the first X-ray telescope called Einstein was put into space thus opening 
the possibility for direct investigation of thermal emission from isolated neutron stars. A signif- 
icant contribution to this study was provided by ROSAT in 1990's. Currently operating observa- 
tories such as Chandra and XMM-Newton have greatly increased the quality and availability of 
observations of thermal radiation from neutron star surfaces. 

The X-ray radiation from an isolated neutron star in general can consist of two distinguishable 
components: the nonthermal emission and the thermal emission. The nonthermal component is 
usually described by a power-law spectral model and attributed to radiation produced in the pulsar 
magnetosphere while the thermal emission can originate either from the entire surface of cooling 
neutron star or the small hot spots around the magnetic poles on stellar surface (polar caps and 
adjacent areas). 

Thermal X-ray emission seems to be a quite common feature of radio pulsars. The black body 
fit allows us to obtain directly the temperature (T s ) of the hot spot. Using the distance (D) to the 
pulsar and the luminosity of thermal emission (Lb i) we can estimate the area (Abb) of the hot spot. 
In the most cases Abb differs from the conventional polar cap area A pc « 6.2 x 10 4 P~ 1 m 2 , where 
P is pulsar period. We use parameter b = A pc /Abb to describe the difference between A pc and Abb- 

1.2 The case with b < 1 

In the most cases observed hot spot area (Abb) is larger then the conventional polar cap area. 
We can distinguish two types of pulsars in this group, with b <C 1 and b <l. 

The first type is associated with observations of thermal emission from the entire stellar surface 
and can be used to test cooling models. Although we have to remember that for young pulsars 
(t = 1 kyr) the nonthermal component dominates, making it impossible to measure accurately the 
thermal flux. As a pulsar becomes older, its nonthermal luminosity decreases faster then the thermal 
luminosity up to the end of the neutrino-cooling era (t ~ 1 Myr). Thus, the thermal radiation from 
the entire stellar surface can dominate at soft X-ray energies for middle-age pulsars (t ~ 100 kyr) 
and some younger pulsars (t ~ 10 kyr) [1]. 

Some observations show that the hot spot area is larger then the conventional polar cap area 
but still significantly less then the area of the star. In this case the radiation comes from the surface 
adjacent to polar cap. Therefore, it can not be explained by cooling of the surface and some 
addition heating mechanism is needed. The model of such heating is based on the assumption 



2 



Partially Screened Gap 



Andrzej Szary 



that the pulsar magnetic field near the stellar surface differs significantly from the pure dipole 
one. The calculations show that it is natural to obtain such geometry of magnetic field lines that 
allows the pair creation in the closed field lines region. The pairs move along closed magnetic field 
lines and heat the surface beyond the polar cap on the opposite side of the star (Fig. 1). In such 
scenario heating energy is generated in IAR (outward particles) hence the luminosity should be 
same order as the nonthermal one but it is hard to predict which of them would prevail in X-ray 
flux of an old neutron star. However, it cannot be ruled out that the thermal one may be dominant 
as suggested by some observations. In most cases area of a such heated surface may be larger (but 
not necessarily) then the conventional polar cap area. That makes the estimation of parameters of 
black-body radiation even harder. More detailed description of this phenomenon can be found in 
[2]. 
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Figure 1: Cartoon of the magnetic field lines in the polar cap region. Red lines are open field lines and 
green dashed lines correspond to the dipole field. The blue arrows show direction of the curvature photons 
emission. 

1.3 The case with b > 1 

In some cases observed hot spot area (Abb) is less then the conventional polar cap area (b > 1). 
The model mentioned in 1.2 can be used in order to explain this phenomenon, because the size of 
modified polar cap is much smaller then size of conventional polar cap (see Fig. 1). The surface 
magnetic filed can be estimated by the magnetic flux conservation law as b = A pc /Abb = B s /Bd. 
Where B d = 2.02 x 10 12 (PP_i 5 ) a5 , P is the pulsar period in seconds and P_i 5 = P/10" 15 is the 
period derivative. 

Medin & Lai [3] calculated the condition for the formation of a vacuum gap above neutron 
star surface. In Fig. 2 we present positions of pulsars with derived surface temperature (T s ) and 
hot spot area (Abb) on the B s — T s diagram where B s is estimated as B s = bBj. Red line represents 
dependence of r cr i t on B s . We can see that in most cases the pulsars positions follow with B s — T cr { t 
theoretical curve. Few cases which do not coincide with the theoretical curve can be explained by 
heating the surface outside the polar cap (see 1.2). 
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Figure 2: Diagram of the surface temperature (T^ = T s /10 6 ) vs. the surface magnetic field (Bu = B S /10 U ). 
The red line is the critical temperature (T^t) evaluated from [3] and dashed lines corresponds to uncertainties 
in calculations. Left panel includes all pulsars with b > 1 while right panel is zoom of the upper right part 
of the graph. Error bars corresponds to la. 



According to our model the actual surface temperature equals to the critical value (T s ~ r cr i t ) 
which leads to the formation of Partially Screened Gap (PSG) above the polar caps of neutron star 
[4]. Hot spot parameters derived from X-ray observations of isolated neutron stars are presented in 
Table 1. 



2. The Model 

The PSG model assumes existence of heavy (Fe 56 ) ions with density near but still below 
corotational charge density (PgjX thus the actual charge density causes partial screening of the 
potential drop just above the polar cap. The degree of shielding can be described by parameter 
r] = 1 — pi/pa where pi is the charge density of heavy ions in the gap. The thermal ejection of 
ions from surface causes partial screening of the acceleration potential drop 

AV = riAV max = 2xriB s h 2 /cP, (2.1) 

where AV max is the potential drop in vacuum gap, h is the gap height, B s = bBj surface magnetic 
field (applicable only if b > 1). Using calculations of Medin & Lai [3] we can express the depen- 
dence of the critical temperature on pulsar parameters as 

r crit = 1.6 x 10 4 (((PP_ 15 )°' 5 &) L1 + 17.7^ (2.2) 

or r crit = 1.1 x 10 6 (b\$ +0.3), where B u = B S /10 U . 

The actual potential drop AV should be thermostatically regulated and there should be estab- 
lished a quasi-equilibrium state, in which heating due to electron/positron bombardment is bal- 
anced by cooling due to thermal radiation (see [4] for more details). The necessary condition for 
formation of this quasi-equilibrium state is 

oT s 4 = T]eAVcn Gh (2.3) 
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where a is the Stefan-Boltzmann constant, e - the electron charge, hqi = Pgj/z= I Ax lO n bP°_^ 5 P~ 05 
is the corotational number density. Using equations (2.2) and (2.3) we can express the acceleration 
potential drop as 

AV = 7.3 x lO 5 ^ -1 ^ 1 (B 1 ^ +0.3)V, (2.4) 
and finally using equations (2.1) and (2.4) we can estimate the gap height in PSG model as 

hr]=6B^(B\i+03) 2 P (2.5) 

As we see both AV and h depend on shielding factor T] which actually depends on the details 
of the avalanche pair production in the gap. First we need to determine which process, Curvature 
Radiation (CR) or Inverse Compton Scattering (ICS), is responsible for the pair production. We 
will need following parameters: Z acc - the distance which a particle should pass to gain Lorentz 
factor equal to y acc , l e - the mean length an electron (or positron) travels before a gamma-photon 
is emitted, and Z p h - the mean free path of gamma-photon before being absorbed by the magnetic 
field. 

2.1 Acceleration path 

In the frame of PSG model, using formalism described in [5], we estimate the component of 
electric field along the magnetic field line in the gap as 

477"/? 

E~n— ^(h-z) (2.6) 
cF 

which vanishes at the top z = h. The Lorentz factor of particles after distance Z acc can be calculated 
as follows 

Y^ = —J Z2 Edz^n A ^{z2-zx){h- Z ^) (2.7) 
m e c l J Z1 m e c 5 P 2 

where m e is mass of an electron and zi — Zi = Zacc- Then we can approximate zi + Z2 ~ 2z and 
assume z~h/2, thus 

m e c 3 P 

'ace — 7acc , • (2.o) 

2ithT]B s e 

2.2 Electron or positron mean free paths 

The mean free path of electron or positron (l e ) can be defined as the mean length that a particle 
passes until a gamma-photon is emitted. In the case of CR electron mean free path can be estimated 
as a distance that particle with Lorentz factor y travels during the time which is necessary to emit 
curvature photon (see [6]) 

P cr V 1 9h q Kc 



le,cr~c[-^ =T-y, (2.9) 

\E YiCr J 4 ye 1 

where P cr = 2y A e 2 c/33i 2 is the power of the CR, E YiCr = 3hy 3 c/23i is the characteristic photon 
energy, 9t - the curvature radius of magnetic field lines. 

For the ICS process, calculation of the electron mean free path l e ^ cs , is not as simple as that 
of the CR process. Although we can define l e j cs in a same way as we defined l 6iCr but it is difficult 
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to estimate a characteristic frequency of emitted photons. We have to take into account photons of 
various frequencies with various incident angles. An estimation of the mean free path of electron 
(or positron) to produce a photon is [7] 



Mi 



Mo ^0 



a f (e,iii)(l-piii)n vh (e)dediii 



(2.10) 



where £ is the incident photon energy in units of m e c 2 , jii = costy is the cosine of the photon 
incident angle, j8 = v/c is the velocity in terms of speed of light, 



4n e 2 



-de 



1 



(2.11) 



" ph A3exp(8/0) 

represents the photon number density distribution of a semi-isotropic blackbody radiation, 9 = 
kT s /m e c 2 , k is the Boltzmann constant, and X c = h/m e c = 2.424 x 10 -10 cm is the electron Comp- 
ton wavelength. Here a' is the cross section of scattering in the particle rest frame. In the Thomson 
regime cross section of scattering in the particle rest frame can be written as [8] 



+ 



(w+1) 2 (u-l) 2 +a 2 _ 



(2.12) 



where Ot is the Thomson cross section, u = £'/£#, a = ^OCf£B, ^ = 7^(1 — Pfr) is the incident 
photon energy in the particle rest frame in units of m e c 2 , £# = B s /B q is the electron cyclotron 
resonance energy in units of m e c 2 , B q = m 2 e c 3 /eh = 4.413 x 10 13 G, and 0Cf = e 2 /he is the fine- 
structure constant. 

Equation (2.12) however does not include the quantum relativistic effects of a strong magnetic 
field (B s > 0.lB q ), therefore, it can be used only at altitudes where magnetic field is relatively 
weaker. For strong magnetic fields we can use an approximation proposed in [9] and then an 
approximate averaged (polarization-independent) cross section can be written as follows 



,_ 3a T ee ,2 {l + ti 2 



2e 



1 



+ ■ 



1 



(e-p q ) 2 (e + A 



(2.13) 



where fi q = B s /B q . Equation (2.13) represents the exact cross section when the particle after scat- 
tering falls on the zero Landau state (see [9] for more accurate results). 

We should expect two modes of ICS process, resonant and non-resonant. The resonant ICS 
takes place if the photon frequency in particle rest frame equals to the cyclotron electron frequency. 
Non-resonant ICS includes all scattering processes of photons with frequencies around the maxi- 
mum of the thermal spectrum. At resonance frequency equation (2.13) has a singularity, so it could 
be used only for non-resonant case. For resonant case, one should use approach proposed in [10], 
thus the cross section in the particle rest frame is 



o' res c^2^--8{e'-e b ) 
m e c 



(2.14) 



Equation (2.14) can be also used for strong magnetic fields (j3 q > 0.1), since the resonant condi- 
tion e' = Eb holds regardless of field strength. Although one have to remember that this equation 
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represents the upper limit for cases with very high magnetic fields (see [9] for more details). The 
particle mean free path above a polar cap for the resonant ICS process is 



levies ^ 



/ / o'nsil-PlidnphWdedHi 
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2xe 2 h 
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dfXi 



(2.15) 



For the altitudes of the same order as the polar cap size we can use jio = 1, fi\ = as incident 
angle limits for outflowing particles, and /lo = 0, ji\ = — 1 as incident angle limits for backflowing 
particles. 

2.3 Photon mean free path 

A photon with energy E r > 2m e c 2 and propagating with nonzero angle \jf with respect to an 
external magnetic field can be absorbed by the field and as a result electron-positron pair is created. 
The optical depth of a such process can be defined as [10] 



(2.16) 



where s p h is a distance traveled by an photon, R^j_ = ± siny/ is the attenuation coefficient for 
the || or _L polarized photons, and R f is the attenuation coefficient in the "perpendicular" frame (the 
frame where the photon propagates perpendicular to the local magnetic field). 

The total attenuation coefficient for pair production is given by R f = Hjk^j & where R'- k is the 
attenuation coefficient for process (channel) in which the photon produces an electron in Landau 
level j and positron in Landau level k, and the sum is taken over all possible states for the electron- 
positron pair. Since pair production is symmetric with respect to the electron and positron, R f k j = 
R' k j for simplicity we will use R f j k to represent the combined probability of creating the pair in either 
(jfc) or (kj) state. For a given channel (jk), the threshold condition for pair production is: 



eL>e',+eL 



(2.17) 



where E' y = Eysini// is the photon energy in the perpendicular frame and E' n = m e c 2 ^/\ + 2j8^n 
is the minimum energy of an electron/positron in Landau Level n. In dimensionless form this 
condition can be written as 



Ey 

2mc 2 



Ey 1 

^- JL 2 Sin ^> o 

2m P c L 2 



\+2P q j+J\+2$ q k 



The first nonzero attenuation coefficients for both polarizations are given in [11]: 
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x > x\ — 1; 



(2.18) 



(2.19) 



^- 2x ££7r=hr**- — ( 1+ \ A ^) /2; (220) 

yX 1 Pq+ 4x 2 



where ao is Bohr radius (let us note that R' ± 00 = and higher orders of attenuation coefficients 
should be used if x > x 3 = (l + ^1+4/3^) /2). 
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The optical depth is defined as: 

T= / R(s)ds= / fl'(s)sini/A<£y (2.21) 
Jo Jo 

We can assume because all high energy photons (x > 1) will produce pairs much earlier 

then \f/ reaches value near unity. In this limit sin \j/ ~ 5ph/9? so relation between x and s p h can be 
expressed by 



9t 2m e c 2 ' 

The equation (2.21) can be rewritten as 



(2.22) 



T = Ti+T|| >2 + Tj_ j 2 + ... 

*i= [ Rioods, %2= I ( R \\w + R \\$\)ds, T ^,2= / R^oids (2.23) 

where si and ^2 are distances which the photon should pass in order to have energy x\ and X2 
respectively (in perpendicular frame of reference). Let us note that 5*1, S2, and S3 are of the same 
order and if s < s\ attenuation coefficient is zero. As shown in [10] for strong magnetic fields 
(Pq ^ 0-1) T i> %2> an d ^±,2 are much larger then one. Therefore, the pair production process 
happens according two scenarios. If fi q > 0.1 photons produce pairs almost immediately upon 
reaching the first threshold, so the created pairs will be in low Landau levels (n < 2). If f5 q < 0.1 
photons will travel longer distances to be absorbed and created pairs will be in higher Landau 
levels. 

Thus, for strong magnetic fields (fi q > 0.1) electron mean free path can be calculated as 

/ P h~*i=*^, (2.24) 

and for relatively weak magnetic fields (j8 g < 0.1) we can use asymptotic approximation derived 
byErber [12] 

I AA h ^ex ( 4 ) (2 25) 

ph (e 2 /hc) m e c f5 q sin \j/ ^ \ 3xf3 q J 

2.4 Partially Screened Gap height 

As we mentioned above, PSG can exist if equation (2.3) is satisfied. On the other hand for the 
heating of the stellar surface the high enough flux of back- streaming particles is required. Thus, 
we need to estimate shielding factor T] and gap height h which are the main parameters of PSG. 
Let us note that T] and h are connected with each other by equation (2.5). Generally gap height can 
by defined as h w l e + l p h with the necessary condition that l e > l acc . Latter corresponds to demand 
that the particle should gain the energy that is required for photon emission by either ICS or CR 
processes. The Fig. 3 shows dependence of free paths on the particles Lorentz factor 7 for some 
particular pulsar parameters (dependence on pulsar parameters will be discussed in section 3). Let 
us note that this free paths do not depend on the gap height h (see equations (2.9), (2.15), (2.24)). 
As it follows from equation (2.8) l acc also does not depend on a gap height because h and T] are 
connected by equation (2.5). 
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Results presented in the Fig. 3 do not allow us to define the gap height unambiguously but we 
can already find that CR is the dominant pair creation process for PSG scenario. Although l e j CS < 
l 6iCr for particles with y ^ 10 3 — 10 4 the gamma-photon production by ICS process is not effective 
because Z acc «l e i CS , i.e. the particles will be accelerated to higher energies (y~ 10 5 - 10 6 ) before 
they would upscatter X-ray photons emitted from the hot polar cap. 




10 3 10 4 10 5 10 3 10 4 10 5 

7 7 



Figure 3: Diagram of the electron/positron mean free path for both CR and ICS processes. The left panel (a) 
corresponds to calculations for B\4 = 1.5, T$ = 1.9, 9?6 = 0.1, P = 1 while the right panel (b) corresponds 
to calculations for Bu = 3.5, Ts = 4.4, 9t 6 = 1, P = 1. Acceleration path was calculated for the same 
gap parameters (h = 20 m, r\ = 0.03) in both panels (see discussion in section 2.4). Although l e j cs < l e ^ cr 
for particles with y ~ 10 3 — 10 4 the y-photon production by ICS process is not effective because for these 
particle Z acc «l e i cs - That means that they will be accelerated to higher energies (y~ 10 5 - 10 6 ) before they 
would upscatter X-ray photons emitted from the stellar surface. 

As soon as we determine the dominant process, i.e. CR, which is responsible for gamma- 
photons emission we can estimate the gap height. Curvature emission by a particle is effective 
for Lorentz factors y ~ 10 5 — 10 6 (when l e ^ cr < Z acc ). The reaction force although is not high 
enough to stop the acceleration by electric field. Equilibrium between acceleration and deceleration 
(by reaction force) would be established when the CR power would equal to "electric power" 
(P cr = P acc , where P acc = A7ir\vB s (h — z)/cP is work done by electric field in unit of time). Using 
the characteristic Lorentz factors of radiating particles we can find characteristic frequencies of 
CR photons and thus, we can check whether the necessary condition (l e + Z p h < h) for cascade 
formation is satisfied. 

The Fig. 4 shows dependence of the CR photon mean free path (l p h), the particle mean free 
path (Zg), and the acceleration path (l aC c) on the energy expressed in units of m e c 2 . Top axis shows 
the Lorentz factors of particles which emit the curvature photons with energy shown on bottom 
axis. While the energy of particle falls in blue region in the Fig. 4 there will be no CR because 
the particle will be accelerated to higher energies before it travels distance enough to emit one 
curvature photon (l e > Z acc ). For the energies in green region CR is most efficient. The particle 
in the PSG never can reach energies falling in red region because reaction force caused by CR 
process is bigger then acceleration force. Thus, the characteristic Lorentz factor y c of particles (and 
characteristic photon energy E c ) corresponds to value defined by the border between green and red 



9 



Partially Screened Gap 



Andrzej Szary 



region. Therefore, we can define gap height for the characteristic energies y c and E c as 



l e (r=Yc) + lph(E r = E c )<h 



(2.26) 



Then we have to calculate l e , l p h for different values of h and check which value satisfies condition 
(2.26). The Fig. 4 shows two cases: panel (a) that corresponds to h = 10 m and panel (b) which 
corresponds to h = 40 m. From panel (a) we can see that condition (2.26) is not satisfied while 
for gap height h = 40 m cascade process will form. This approach allows us to calculate minimum 
height at which gap sparking breakdown is possible and thus, we can estimate unambiguously the 
gap height in the frame of the PSG model. 
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Figure 4: Diagram of the photon mean free path, the particle mean free path, and the acceleration path vs. 
the photon energy and the Lorentz factor of accelerated particles (top axis) for specific pulsar (Bu = 3.5, 
Ts = 4.39, = 0-3, P = 1). See description of blue, green and red regions in the text above equation (2.26). 
The black horizontal line corresponds to the gap height. Calculations for panel (a) were done using h = 10 
m, 77 = 0.03 while for panel (b) using h = 40 m, 77 = 0.01. 



3. Observational consequences 

In order to compare the model with observations we need to define main parameters of PSG 
model for a given pulsar with known period and period derivative. 

3.1 PSG model parameters 

We can distinguish two types of PSG parameters: observed and derived. As we mention above 
in some cases when X-ray observations are available we can directly estimate surface magnetic 
filed B s . On one hand B s can be calculated using the size of the hot spot Abb, and on the other hand 
we can find B s using estimation of the critical temperature and assumption that T s = T cr [ t . One of 
the most important requirements of PSG model is that these two estimations should coincide with 
each others. As it is clear from Fig. 2 in most cases when the hot spot parameters are available 
this requirement is fulfilled. Thus, we can assume the characteristic values of B s vary in range of 
(1 — 4) x 10 14 which corresponds to critical/surface temperature in range of (1.1-4.5) x 10 6 (see 
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Table 1). Using these values we can estimate derived parameters of PSG such us the gap height 
h, the shielding factor T] and the Lorentz factor of primary particles y c . Let us note that these 
parameters depend also on the curvature radius of magnetic field lines 9?. The curvature can not 
be neither observed or derived but modeling of surface magnetic field (see Fig. 1) indicates that 
the curvature radius varies in the range of (0.1 — 1) x 10 6 cm. Below we will discuss influence of 
pulsar parameters such as the magnetic field, the curvature of field lines and the period on derived 
PSG parameters. 

3.2 Influence of the magnetic field 

The conditions in PSG are mainly defined by the surface magnetic field. In Fig. 5 panel (a) 
we present dependence of the gap height on the surface magnetic field calculated according to 
approach described in section 2.4. It is clear that the gap height decreases as the surface magnetic 
field increases. The Fig. 5 panel (b) shows dependence of the shielding factor on the surface 
magnetic field calculated using equation (2.5). We can see that for stronger magnetic fields T] 
increases which means that the density of heavy ions above the polar cap decreases. Let us note 
that the surface temperature (T s ) stays very near to the critical temperature (T cr { t ) which is shown 
on the top axis of the diagrams. In Fig. 5 panel (c) we present dependence of Lorenz factors of the 
particles accelerated in PSG. The Green line (70) presents the values corresponding to the boundary 
between blue and green region in the Fig. 4 while the red line corresponds to the boundary between 
the green and the red region which is the characteristic value (y c ) for primary particles. We see that 
y c is very slightly affected by the magnetic field strength. 

One can expect that for higher temperatures (which correspond to stronger magnetic fields) 
the gap breakdown could be dominated by ICS process. However this is not the case because the 
particle mean free path is much higher then the acceleration path (l e j cs » Z acc ) even for strong 
magnetic fields. 

T 6 T 6 T 6 

23456 23456 23456 




Figure 5: Dependence of the gap height (panel a), the shielding factor (panel b), and the particles Lorentz 
factor (panel c) on the surface magnetic field. Solid lines correspond to calculations for 9^ = 0.1 while 
dashed lines correspond to calculations for 9t 6 = 1 . For both cases the pulsar period P = I was used. The 
green line presents the values corresponding to the boundary between blue and green region in the Fig. 
4 while the red line corresponds to the boundary between green and red region. Corresponding critical 
temperature is shown on top axis of diagrams. 
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3.3 Influence of the field lines curvature radius 

The curvature of magnetic field lines significantly affects the PSG parameters since the CR 
process is responsible for PSG breakdown. As it can be expected the gap height decreases for 
smaller curvature radius of magnetic field lines (Fig. 6 panel a) since CR is more efficient for 
smaller values of the curvature radius. Consequently the shielding factor decreases with increasing 
curvature radius of magnetic field lines (Fig. 6 panel b). 

The Lorentz factor of primary particles increases for bigger radius of curvature (Fig. 6 panel 
c), since both conditions l e ~ Z acc and P cr ~ P acc are satisfied for higher Lorentz factors. 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

3^6 3^6 3^6 

Figure 6: Dependence of the gap height (panel a), the shielding factor (panel b), and the particles Lorentz 
factor (panel c) on the curvature radius of magnetic field lines. Solid lines correspond to calculations for 
Bu = 3.5 while dashed lines correspond to calculations for Bu = 1.5. For both cases the pulsar period 
P = 1 was used. 



3.4 Influence of the pulsar period 

As we can see from figure 7 panel (a) and panel (c) neither the gap height nor the Lorentz 
factor of primary particles depend on pulsar period. This is one of the most important result of our 
calculations. The shielding factor is just proportional to period (Fig. 7 panel b) as is expected from 
equation (2.5). 



1 1 1 1 1 — I I — l 1 1 1 1 — I F — I 1 1 I I ~ [ 




P[s] P[s] P[s] 



Figure 7: Dependence of the gap height (panel a), the shielding factor (panel b), and the particles Lorentz 
factor (panel c) on the pulsar period. Solid lines correspond to calculations for B\4 = 3.5 while dashed lines 
correspond to calculations for Bu = 1.5. For both cases the radius of curvature 3i^ = 0.l was used. 
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4. Conclusions 

X-ray observations of pulsars show that temperature of the hot spot is about few million kelvins 
while its area is much smaller then the conventional polar cap surface which is calculated assuming 
purely dipolar geometry of pulsar magnetic field. Such observations can be easily explained in the 
frame of PSG model which assumes that the geometry of magnetic field near the stellar surface 
differs significantly from dipolar one, actually the field is much stronger and curved. At the same 
time the surface temperature is about the critical temperature which depends only on the magnetic 
field strength. 

In this paper we calculated dependence of the gap parameters such as the gap height, the 
shielding factor, and characteristic energies of primary particles on the surface magnetic field, the 
curvature radius of field lines, and the pulsar period in the frame of PSG model. The surface 
magnetic field can be calculated using the X-ray observations (if such data are available), however 
the curvature radius can be only estimated by simulations of different geometry of magnetic field 
lines. Let us note that the gap height is the most important parameter for all models of the Inner 
Acceleration Region in pulsars. We can define the gap height as a sum of the particle mean free 
path and the photon mean free path, l e + l p h- In order to estimate the PSG height we discuss 
two processes responsible for cascade pair production: the Inverse Compton Scattering and the 
Curvature Radiation. We found that CR is the dominant process for the sparking breakdown of the 
PSG. Although ICS is more efficient for particles with y ^ 10 3 — 10 4 , the particles are accelerated 
to higher energies 7 ~ 10 5 — 10 6 before they upscatter X-ray photons emitted from the polar cap. 
As soon as particles Lorentz factor 7 > 10 5 CR is more efficient then ICS. In section 3 we show 
the dependence of PSG parameters (h, 77, y c ) on the pulsar parameters (B s , 91, P). Since CR is the 
dominant process for gamma-photon emission the PSG parameters strongly depend on curvature 
radius of magnetic field lines. Pulsars with smaller curvature radius should have the lower gap 
height and also the higher shielding factor, consequently the density of heavy ions should be lower. 
In pulsars with stronger surface magnetic field gap heights should be smaller and also shielding 
factor should be smaller. Our calculations show that on one hand h and y c do not depend on pulsar 
period, but on the other hand T] increases along with increase pulsar period. The evaluated PSG 
parameters should play decisive role in pulsar emission models. 

We found that PSG model is suitable to explain both cases: when the hot spot is smaller then 
the conventional polar cap, and vice versa when the hot spot is larger then conventional polar cap. 
In the latter case the surface is heated by particles created in the closed magnetic filed lines region 
by photons emitted in the open magnetic filed lines region. Let us note that in the purely dipolar 
magnetic field geometry photons emitted tangent to magnetic field lines always stay in the open 
field lines region. 
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